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1. INTRODUCTION 

Consider a geodesic line between two points A and B on 
the surface of the Earth. Given the position of A, the length of 
the line and its azimuth at A, we wish to determine the posi- 
tion of B and the azimuth of the line there. This problem oc- 
curs so frequently that I undertook to construct tables to sim- 
plify the computation. In order to explain the method clearly, I 
start by deriving the fundamental properties of geodesic lines 
on a spheroid of revolution. Even though aspects of this 
derivation may already be well known, the benefit of having 
the entire development presented together outweighs the cost 
of repeating it. 1 



2. THE CHARACTERISTIC EQUATION FOR A GEODESIC 

Take two points A and B on the surface on a spheroid 2 
of revolution joined by some specified curve. Consider two 
neighboring points on the curve with latitudes <f> and <f> + d<j> 
and longitudes relative to A of w and w + dw (measuring east 
positive). Let the distance between them be ds, the azimuth of 
line directed toward A be a (measured clockwise from north), 



the radius of the circle of latitude be r, and the meridional 
radius of curvature by R; then we find 3 



"This is an English translation of (jber die Berechnung der geographischen 
Langen und Breiten aus geoddtischen Vermessungen. Astronomische Nach- 
richten 4(86), 2 41-254 (1826), doi[10 .1002/asna.l826004160 i] Th e paper 
also appears in Abhandhmgen von Friedrich WUhelm Bessel, Vol. 3 pp. 5— 
14 (W. Engelmann, Leipzig, 1876). The translation has been prepared and 
edited by Charles F. F. Karney i ckarney@sarnoff.com > and Rodney E. Deakin 
irod.deakin@rmit.edu.au i, with the assistance of Max Hunter and Stephan 
Brunner. The mathematical notation has been updated to conform to cur- 
rent conventions and, in a few places, the equations have been rearranged 
for clarity. Several errors have been corrected, a figure has been included, 
and the tables have been recomputed. A transcription of the original paper 
with the updated mathematical notation and with the corrections is available 
at arXiv:0908.1823 A contemporary, but partial, translation into English ap- 
peared in Quart. Jour. Roy. Inst. 21(41), 138-152 (1826). 

1 In Secs.ffjg] Bessel gives a concise summary of the work of several other 
authors, notably, Clairaut, Euler, du Sejour, Legendre, and Oriani. Bessel's 
contributions, which start in Sec. [5] consist of his methods for expanding 
the distance and longitude integrals and his compilation of tables to provide 
a practical method for computing geodesies. Two sentences have been 
omitted from this translation of the introduction. In one, Bessel refers to 
two letters he published earlier in the Astronomische Nachrichten which 
do not, however, have a direct bearing on the present work. In the other, 
he criticizes "du Sejour's method," but without providing details; in any 
case, such criticism is misplaced because du Sejour had died over 30 years 
earlier and Bessel does not cite more recent work. 

2 "Spheroid" here is used in the sense of a shape approximating a sphere. 
Sections f5]and[3]treat the case of a rotationally symmetric earth. In Sec. [4] 
Bessel specializes to a rotationally symmetric ellipsoid. 



dv 

cos ads = —Rdcj) = — — -, 
sin <p 

sin ads = —rdw, 



(D 



which gives 



ds = y/R 2 d(j) 2 +r 2 dw 2 . 



If we write p for dqb/dw and U for y/ R 2 p 2 + r 2 , this becomes 
ds = U dw. 

The distance along the curve between the two points A and B 
is therefore 

Udw, 

where the integration is from A to B. If the curve is the 
geodesic or shortest path, then the relation between <p and w 
must be such that the integral is a minimum. If we perturb this 
relation so that <p is replaced by <f> + z where z is an arbitrary 
function of w which vanishes at the end points (because these 
points lie on both curves), then the perturbed length, 

s = I U' dw, 



must be larger than s for all z. 

Expanding U(cf>,p) in a Taylor series, we obtain 4 

U' - U +—Z+-— ^ 

d(j) dp dw 

and therefore we have 

dU z+ dU dz \ dw + 
def) dp dw J 

where we have explicitly included terms only up to first order 
in z. For s to be a minimum, we require that 

dU_ z dU dz_\ w+ >Q 
d(j> dp dw J 



3 The minus signs appear in {T} because a is the back azimuth, pointing to 
A, while ds advances the geodesic away from A. In this section, Bessel 
assumes an easterly geodesic so that ds /dw > 0. However the final result, 
Eq. {2)> is general. 

4 The notation here employs partial derivatives instead of Bessel's less for- 
mal use of differentials. 
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for all z. Since this must also hold if z is replaced by — z 
and since we can take z so small that the first order terms are 
bigger that the sum of all the higher order terms (except if the 
first order terms vanish), it follows that the condition that s be 
minimum is 

(dU_ 9U^±\ dw _ {) 
\dcj) dp dw J 

Integrating the second term by parts to give z(dll /dp) — 
J z [d(dU /dp) /dw] dw and remembering that z vanishes at 
the end points, we obtain 



r au 
z \dj 



_d_fdU_ 

dw \ dp 



dw = 0. 



Since this integral must vanish for arbitrary z, we find 5 

&U_ _d_fd£\ (| 

d<j) dw \ dp 



d fdU 



or, multiplying by dip / dw = p, 

dU # dU dp dp dU 
dcf) dw dp dw dw dp ^ dw \ dp 

which on integrating with respect to w becomes 



IT ( dU 



const. 



Substituting r 2 + R 2 p 2 for U, we obtain 6 

r 



y/1 + {R 2 /r 2 )p< 



= — rsina = const., 



which is the well known characteristic equation of the geo- 
desic. 

If the azimuth of the geodesic at A (in the direction of B) is 
a' and the distance of A from the rotation axis is r', we have 

r sin(a' + 180°) = rsina, 




Figure 1 Spherical triangles on the auxiliary sphere. EAB is the 
geodesic, N is the pole; EFG is the equator; and NE, NAF, and 
NBG are meridians. 



and equation (f2j) becomes 



cos it sm a 



cos usm a. 



(3) 



This equation relates two sides of a spherical triangle, 8 90° — 
v! and 90° — u, and their opposite angles, 360° — a and a'. 
The third side a and its opposite angle w will appear in the 
following calculations giving elegant expressions for the joint 
variations of s, u and w. In particular, using the well known 
differential formulas of spherical trigonometry, we find 9 



du = — cos a da, 
cos udui = — sin a da. 



Substituting these in equations (HJ and expressin 
of u gives 



-rsma. 



(2) 



smu 
ds = a— — -da, 



sin it 
dw = — — -dw. 



lg r in terms 



(4) 



3. THE AUXILIARY SPHERE 



4. THE EQUATIONS FOR A GEODESIC ON AN ELLIPSOID 



Let the maximum distance of the spheroid to the rotation 
axis be a, so that r and r' are less than or equal to a; we can 
then write 7 



I now assume that the meridian is an ellipse with equa- 
to rial rad ius a, polar semi-axis b, and eccentricity e = 
Va 2 — b 2 /a. w The equation for an ellipse expressed in terms 



a cos it 



a cos it, 



5 This is the Euler-Lagrange equation of the calculus of variations. 

6 A. C. Clairaut gives a geometric derivation of this result in Mem. de 
l'Acad. Roy. des Sciences de Paris 1733, 406^116 (1735) The equation 
also follows from conservation of angular momentum for a mass sliding 
without friction on a spheroid of revolution. 

7 The quantity u is the reduced or parametric latitude. 



8 See the triangle ABN on the "auxiliary sphere" in Fig.fT] Equation (5) is 
the sine rule applied to angles A and B of the triangle. 

9 Here and in the rest of the paper, the differentials give the movement of 
point B along the geodesic denned with point A and a' held fixed. 

10 In Bessel's time, it was known that the earth could be approximated by an 
oblate ellipsoid, a > b, but the eccentricity had not been determined accu- 
rately. Therefore, Bessel computes tables which are applicable to ellipsoids 
with a range of eccenU'icities. 
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of cartesian coordinates is 

2 2 

_ +yi = 1 

a 2 6 2 

Differentiating this and setting dy/dx — — cot </>, we obtain 
xsin0 ?/cos< 



6 2 



= 0; 



eliminating y between these equations then gives 

a cos <$> 



yl - e 2 sin 2 </> 

The quantity a; is the same as r = a cos u, which gives the 
relationships between <fi and u, 



cos< 



cos u 



e 2 sin 2 <A 



COS( 



sin <fiy/l — e 2 



e 2 sin 2 d> 



tan u = tan §\/ 1 — e 2 , tan = 



cos it Vl — e 2 
\/l — e 2 cos 2 u 

sin it 
\/l — e 2 cos 2 u 
tanu 



and 



sin it 
sin rf> 



e 2 cos 2 u. 



Substituting this into (01, we obtain the differential equations 
for a geodesic on an ellipsoid 



ds = ay 1 — e 2 cos 2 u da, 



dw = \/\ — e 2 cos 2 u du). 



(5) 



5. THE DISTANCE INTEGRAL 

To integrate the first of these differential equations, I use 
the three relations between u', u, a', a and cr, 11 

sin u = sin u cos a + cos u cos a 1 sin a, 

— cos u cos a — — sin u' sin a + cos u' cos a cos cr, (6) 

— cos u sin a = cos v! sin a' . 

It is convenient to write these in terms of the auxiliary angles 
m and M defined by 12 

sin u = cos m sin M, 
cos v! cos ol = cos m cos Af, (7) 
cos v! sin a' = sin m. 



11 Referring to Fig.[TJ consider two central cartesian coordinate systems with 
the xy plane containing the geodesic EAB, and either A or B lying on 
the x axis. Equations (6) give the transformation between the coordi- 
nates of N in the two systems, [sin u' , cos u' cos a', cos u' sin a'] and 
[sin u, — cos u cos a, — cos u sin a] , namely a rotation by a about the z 
axis. 

12 The auxiliary angles m and M are an angle and a side of the spherical 
triangle EAN shown in Fig.[T] Equations jTJ are the sine rule on angles E 
and F of triangle AEF, the cosine rule on angle F of triangle AEF, and 
the sine rule on angles A and E of triangle ANE. 



Equations (O then become 13 

sin u = cos to sin( M + cr) , 
cos it cos a = — cos77icos(Af + a), (8) 
cos m sin a = — sin m. 

This gives 

cos 2 tt = 1 — cos 2 m sin 2 (M + cr) , 
and the equation for ds becomes 



ds = aVT - e 2 Jl + k 2 sin 2 (M + cr) At, (9) 



where 



fc = 



e COS 777 



This differential equation may be integrated in terms of the 
elliptic integrals introduced by Legendre. 14 Because the tools 
to compute these special functions are not yet sufficiently ver- 
satile, we instead develop a series solution which converges 
rapidly because e 2 is so small. We readily achieve this by de- 
composing the term under the square root into two complex 
factors, namely 15 



ds = 



VT^ 2 

! 

1 - e 



dax 



yjl - e exp(2i(A7 + a))^l-e exp(-2i(M + cr)) 



where 



_ Vi + fc 2 - l _ 2^1 

6 — — ~~ , rZ 



VT+k 2 + 1 



1 - e 



Expanding the two factors in the radicals in infinite series and 
multiplying the results gives 16 



■J\ _ e 2 

ds = a- da \A - 2B cos 2(M + a) 

1 — e 

- 2Ccos4(A/ + cr) - 2Dcos6(M + a) - ...], 



13 These are analogs of Eqs. jTJ with meridian replaced by NBG. 

14 A. M. Legendre, Exercices du calcul integral Vol. 1 (Courcier, 1811). 

15 The notation has been simplified here compared to Bessel's original for- 
mulation in which k and e are expressed in terms of E through k = tan E 
and e = tan 2 \E. By using t as the expansion parameter and by dividing 
out the factor 1 — e, Bessel has ensured that the terms that he is expanding 
are invariant under the transformation e — > —t,M + u — > n/2 — (M+cr). 
This symmetry causes half the terms in the expansions in e to vanish. 

16 The use of complex exponentials facilitates the series expansions by avoid- 
ing the need to employ awkward trigonometric identities. If we write 

V 1 2 2-4 2-4-6 2-4-6-8 "' Z J j"j J - > 

then the coefficient of cos (2l(M + cr)) e'+ 2 J is a 2 for Z = and 2aj aj + ; 
for Z > 0. 
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where A, B, C, . . . are given by 

n 1 1-1 1 3 1-1-3 1-1 5 

B = -e e A e 5 

2 2-4 2 2-4-6 2-4 

1-1-3-5 1-1-3 7 

~ 2-4-6-8 2-4-6 £ ~ '"' 

„ 1-1 2 1-1-3 1 4 1-1-3-5 1-1 fi 

C = e 2 e 4 e 6 

2-4 2-4-6 2 2-4-6-8 2-4 

M-3-5-7 1-1-3 8 

~ 2-4-6-8-10 2-4-6 € 
_ 1-1-3 3 1-1-3-5 1 5 M-3-5-7 1-1 7 

~~ 2-4-6 £ ~ 2-4-6-8 2~ £ ~ 2-4-6-8-10 2~4 £ 
1-1-3-5-7-9 1-1-3 9 

~ 2-4-6-8-10-12 2-4-6 £ 

etc. 

Integrating the equation for ds starting at a = 0, we obtain 

s = — ^— \Ao - jB cos(2M + a) sin a 
1 — e 

- |C cos(4M + 2<t) sin 2er 

- §L> cos(6M + 3cr) sin 3ct 
-...]. (10) 

6. SOLVING THE DISTANCE EQUATION 

The series ( TTOb gives the distance s between A and B in 
terms of u', a', and tr; if, however, s and a' have been 
measured and u' is known from the latitude at A, then a 
is obtained by solving (TTOb . The latitude of B and the az- 
imuth of the geodesic there are found from (0. Equation (TTOb 
can be solved either by reverting the series or by successive 
approximation — the latter way is however the simplest if the 
tables I have compiled are used. 

I write 17 

<t = — s+ i 8cos(2M+cr)sinCT+7cos(4M+2cr)sin2<7 

+ S cos(6M + 3cr) sin 3ct + . . . , (11) 

where 



648 000 1 - e 



a = 


TT A 


P = 


648 000 2B 


TT A ' 


7 = 


648 000 C 


TT A' 


5 = 


648 000 2D 


TT 3 A' 



etc. 



The units for a, a, 0, ... are arc seconds. Bessel here adopts a conflict- 
ing notation for the coefficient a which should not be confused with the 
azimuth. 



The tables give the logarithms 18 of a, (3, and 7 as a function 
of the argument 

e cos m 

logfc = log 

V 1 — e 

By this choice, the variation of log f3 and log 7 are very close 
to two and four times that of the argument, which simplifies 
interpolation into the table. 19 

We take as/b as the first approximation of a, substitute this 
into the second term to obtain a second approximation, with 
which we recalculate the second term and add the third. The 
convergence of the series is sufficiently fast that, even if the 
argument is 1.1 (which is only possible if the flattening of 
the ellipsoid, 1 — b/a, exceeds jig), the approximation never 
needs to be carried further in order to keep the errors in a 
under 0.001". The term involving 6 does not exceed 0.0005" 
at this value of the argument. 

7. ACCURACY OF THE TABLES 

The values of log a in the table are given to 8 decimal 
places. An error of half a unit in the last place results in an 
error of only 0.0005" or 0.008 toise over a distance corre- 
sponding to a = 12° '4' or 700000 toises. 20 Similarly, I retain 
only sufficient digits in the tabulation of log j3 to ensure that 
the error in this term is less than 0.0005"; for this purpose, I 
use 6 digits at the end of the table and fewer digits for smaller 
values of the argument. The third term never exceeds 0.17", 
even at the end of the table; therefore I include only 3 decimal 
places for log 7. Thus the errors are 0.001" for distances up to 
700000 toises; even if the distance is of the order of a quarter 
meridian (i.e., a — 90°), the error is less than 0.01". 



8. AN EXAMPLE 

In order to illustrate the use of the tables, I consider the 
results from the great survey by von Muffling. 21 Relative to 



18 In this paper, logx denotes the common logarithm (base 10) and we use 
colog x = log(l/x). The tables in the original paper contained a number 
of errors of one unit in the last place. These errors do not, for the most 
part, affect the results obtained from the tables when rounded to 0.001". 
In addition, there were systematic errors in the tabulated values of log 
equivalent to a relative error of order t 2 in which result in discrepancies 
from 1 to 17 units in the last place on the final page (the 6-figure portion) 
of the tables. In calculations involving logarithms, a bar over a numeral 
indicates that that numeral should be negated, e.g., log 0.02 fa 2.3 = 
(—2) + 0.3. In the original paper, logarithms are written modulo 10, e.g., 
log 0.02 fa 8.3. The notation "(— )" in these calculations indicates that the 
quantity whose logarithm is being taken is negative. 

" The columns headed A give the first differences of the immediately pre- 
ceding columns and aid in interpolating the data. Bessel would have used 
a table of "proportional parts" to compute the interpolated values. 

20 The toise was a French unit of length. It can be converted to meters by 
ltoise = 8641igne, 443.296 ligne = lm, or ltoise « 1.949m. 

21 F. C. F. von Muffling, [Asiron. Nachr. 2(27), 33-38 (1824)] 
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Seeberg (point A), the distance and azimuth to Dunkirk (point 

B) are 22 

logs = 5.478 303 14, 
a 1 = 274° 21' 3.18". 

I assume the latitude of the Observatory at Seeberg to be 
4>' = 50° 56' 6.7" and the ellipsoid parameters to be log b = 
6.513 354 64, log e = 2.90 5 4355. 23 
From tanu' = \/l — e 2 tan^', we find 

logtan</>' = 0.090 626 65 



log \7l-e 2 = 1.998 590 60 



log tan u' = 0.089 217 25; u = 50° 50' 39.057". 

Given u' and a', we can compute M, cosm and sinm from 
equations ©I 24 

log sin u' = 1.889 543 51 
log cos v! = 1.800 326 27 
log cos a' = 2.880 037 33 
log sin a' = 1.998 746 62(-) 



log(cos to sin Af) = 1.889 543 51 
log(cos m cos M) = 2.680 363 60 

log sinm = 1.799 072 89(-) 



M = 86° 27' 53.949"; 2M = 172° 55' 47.9" 
log cos m = 1.890 370 63 4M = 345° 51' 36". 

The argument in the tables, log((e/Vl — e 2 ) costo), is 



log 



= 2.906 845 



log cos to = 1.890 371 



Argument = 2.797 216. 
Looking up log a in the tables, and calculating as/b gives 25 
log a = 5.313 998 92 
colog 6 = 7.486 645 36 
logs = 5.478 303 14 

a 



log — = 4.278 94742; 
b 



-s = 5° 16' 48.481". 



22 Seeberg: 50°56'N 10°44'E; Dunkirk: 51°2'N 2°23'E. 

23 In present-day units, this is a 6377 km, flattening / 1/308.6, s ~ 
586 km. In this example, Bessel uses the toise as the unit of length and the 
second as the unit of arc. 

24 Bessel solves 3 equations for 2 unknowns M and m. The redundancy 
serves as a check for the hand calculation and can also improve the accu- 
racy of the calculation, for example, in the case where sin m PS 1. 

25 It is necessary to use second differences when interpolating in the table for 
log a. The argument, 2.797 216, lies q = 0.7216 of the way between 
2.79 and 2.80. Bessel's central 2nd-order interpolation formula for the last 
6 digits of log a gives 401 284 + g(-1941) + \q(q - 1)(1853 - 1004 - 
1028) = 399 892. For the other table look-ups, linear interpolation using 
first differences suffices. 



Adopting this as the first approximation to the value of a, we 
obtain the second by adding the first term in the series ( fTTT i, 

log ,3 = 2.305 94 
logcos(2M + o-) = 1.999 79(-) 
log sin g = 2.963 91 



1.269 64(-) = -18.61". 

We now update the value of this term with the second approx- 
imation of a = 5° 16' 48.5" - 18.6" = 5° 16' 29.9" and so 
obtain as the third approximation: 

log/3 = 2.305 94 
logcos(2M + cr) = 1.999 79(-) 
log sin cr = 2.963 48 



1.269 21(-) = -18.587", 



log 7 = 2.394 
logcos(4M + 2<r) = 1.999 
log sin 2cr = 1.263 



3.656 = +0.005". 

Gathering the terms in ([TTJ gives a = 5° 16' 48. 481" - 
18.587" + 0.005" = 5° 16' 29.899" and so, finally, we de- 
termine a, u and <f> from equations ((H), 

M + ct = 91° 44' 23.848" 
logsin(M + cr) = 1.999 799 71 
log(- cos(A/ + ct)) = 2.482 349 32 
log cos to = 1.890 370 63 
log(-sinm) = 1.799 072 89 



log sinu = 1.890170 34 
log(cosMCOsa) = 2.372 719 95 
log(cosusina) = 1.799 072 89 



log cot a = 2.573 647 06; a = 87° 51' 15.523" 
log cos u = 1.799 377 50 
log tan u = 0.090 792 84 



colog \f\ 



0.001409 40 



log tan cf) = 0.092 202 24; 



51° 2' 12.719'' 



In this example, I carried out the trigonometric calculations to 
8 decimals; however the tables of log a, log/3, and log 7 in 
fact allow a and <j> to be determined slightly more accurately 
than this. If only standard 7 -figure logarithm tables are avail- 
able, the last digits in the tabulated values of log a, log j3, and 
log 7 may be neglected. 
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9. THE LONGITUDE INTEGRAL 

We turn now to the determination of the longitude differ- 
ence w by integrating (|5), 



dw — yl — e 2 cos 2 uduj. 

This integral contains two separate constants m and e, which 
cannot be combined. Thus it not possible to construct tables to 
allow a rigorous solution of this problem which are valid for 
arbitrary e. 26 However, we can achieve our goal by sacrificing 
strict rigor and by making an approximation which results in 
errors which are inconsequential in our application. 
We start by writing 



dw = du> — (l — \/ \ — e 2 cos 2 ujdui, 



and substitute in the second term 



dw 



■ da. 



On integrating, we obtain 



w = oj — sin a cos u 



, /"i-vi 



e'- cos^ u 



■ da. 



Let us write 



1 - VT - e 2 cos 2 u e 2 , , .„._, . 

a = -77 1 + e 2 pcos 2 u)»(l + y); 

cos z u I 

in other words, we set 



l + y = 



2(1 - Vl - e 2 cos 2 u) 
e 2 cos 2 u(l + e 2 p cos 2 u) q 



1 + |;e 2 cos 2 u + |e 4 cos 4 u + ^e 6 cos 6 u + . . . 

(1 + qpe 2 cos 2 u + ^j^- jpe 4 cos 4 u \ 
+ q( % 1 l {q 3 - 2) p 3 e 6 cos 6 U + ..J 

The first three terms in the denominator and in the numerator 
are equal, provided that 

P = -|, ?=-s, 

which gives 

1 + \ e 2 cos 2 -u + |e 4 cos 4 u + ^e 6 cos 6 u + . . . 



1 + y 



1 + |e 2 cos 2 m + |e 4 cos 4 it + ^e 6 cos 6 u + 



1 + -^e 6 cos 6 u + 



26 As a practical matter, it would have been impossible for Bessel to provide 
a complete tabulation of a function of two parameters. He could have tab- 
ulated the function for a fixed value of e, which would greatly reduced 
the utility of his method, especially given the uncertainties in the measure- 
ments of e. Instead, Bessel manipulates the expression for dw to move 
the dependence on the second parameter into a small term that may be 
neglected. 



From this, we see that neglecting y results in an error of or- 
der e 8 or an error in w of ^eV. This would not be dis- 
cernible even in the calculation of long geodesies to 10 deci- 
mal places. 27 

Thus, for the present purposes, we may take y sw enabling 
us to tabulate the integral in a way that is valid for all e. 



10. SERIES EXPANSION FOR LONGITUDE 

Introducing this approximation, we have 

= 2 r da 



w ~ uj — sin m 

1 — | e 2 cos 2 u 



oj — sin m 



da 



If we set 



1 - |e 2 + |e 2 cos 2 msin 2 (M + a) 
, yf\ e cos m 



0" 



3 P 2 



we can express the integral in the second term as 

da 



^/l-|e 2 ^/l + fc' 2 sin 2 (M + <7) 

Following the same procedure used in expanding the integral 
for ds in Sec. [5] we introduce e' defined by 28 



, _ VTTfc 72 -! , _ 2V7 
VTTk^ + i l-e' 

and separate the integrand into two complex factors, 



m - e' exp(2i(M + a)) - e' exp(-2i(Af + a)) 
If we expand these in infinite series, the product becomes 29 
2 



3p2 



(a 1 + ft cos 2(M + a) + 2i cos 4(M + a) 
+ 3(5' cos6(M + <t) + ...W 



27 For a flattening of , the error in the longitude difference over a distance 
equivalent to a quarter meridian, i.e., 10 000 km, is less than 0.000 05". 

28 Bessel gives the relationship between k' and e' in terms of E', where k' = 
tan E' and e' = tan 2 \E' . 

29 There are a series of errors in the original paper leading up to 1121 . Here 
we assume that the original Eq. \ defines a' , 7', . . . , which makes 
this equation analogous to (TT) , and correct the preceding equations to be 
consistent. 
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where 30 



i 



3 1 6 



1 , 1-4 1 , 3 

-e H e 6 

3 3-6 3 

1^4,2 1-4-7 1 /4 

3-6 £ 3-6-9 3 £ 



3-6 y 

1-4-7 1-4 



3-6-9 3-6 



1-4-7 



1-4-7-10 1-4 
' 3-6-9-12 jTtj' 
1-4-7-10 1 



3-6-9 



3-6-9-12 3 



1-4-7-10-13 1-4 , 7 
3-6-9-12-15 3^6 £ 



etc. 



Integrating from cj = then gives 



e sin m 



U> W CJ - 



0" 



+ /3' cos(2M + cr) sin a 

-7'cos(4M + 2cr)sin2cr 

-5'cos(6M + 3(r)sin3c7 + ...V (12) 



11. COMPUTING THE LONGITUDE DIFFERENCE 

The first two coefficients of this series are given in the 4th 
and 5th columns of the tables 31 as functions of the argument 



log k' = log 



: cosm 



The convergence is commensurate with the 3 first columns 
of the tables. We calculate uj using one of the formulas for 
spherical triangles (Sec. [3]l, either 32 

sin a sin a' — sin a sin a sin a sin m 



cosu 



cos u' cos u cos u' 



or 



33 



, sinifu — u) , . 

tan tjCj = f- cot ^(a + a) 

2 cosi(u' + u) 2 

cosi(u'-u) , 

= f- cot t (a — a). 

sini(u' + u) 2V ' 



30 c 



1 See footnotefTSland set (1 - a;)" 1 / 3 = 1 + \x + ±^x 2 + j^x 3 + 
1-4-7-10 4 i 
3-6'9-12' i ' 

The value of /?' in the tables includes the factor of 648 000/-7T necessary to 
convert from radians to arc seconds. 

32 The first two relations are the sine rule for angle N of triangle ABN of 
Fig.fJJ The last relation is obtained, for example, by substituting for sin a' 
from (7). 

33 These are Napier's analogies for angle N of triangle ABN. 



and evaluate w by means of the tables. 

I will continue with the example in Sec.[8]and calculate the 
longitude difference between Dunkirk and Seeberg using this 
prescription. Solving the spherical triangle for lu gives 

log sin a = 2.963 483 83 
log(-sina) = 1.999 695 39(-) 
cologcosu' = 0.199 673 73 



logsin uj = 1.162 852 95(-); u = -8° 21' 57.741". 
The argument for the last two columns of the tables is 

1 °g((v / l e /v / l 3 le 2 ) cosm), giving 
Vie 



log 



2.844 022 



log cosm = 1.890 371 
Argument = 2.734 393. 
Computing the terms in the series dTZb gives 

log a' = 1.698 758 
log(-sinm) = 1.799 073 



log 



3p2 



= 3.811575 



lag rj = 4.278 523 



1.587 929 = +38.719", 



and 



log/3' = 1.703 
log(— sinm) = 1.799 



log 



3 P 2 



= 3.812 



log(cos(2M + a)sina) = 2.963(— ) 

2.277(-) = -0.019". 

The sum of both terms is +38.700", and adding this to lu, we 
find the longitude difference, 

w = -8° 21' 19.041". 

12. CONCLUSION 

This illustration of the use of these tables shows that the 
accuracy of the calculation is limited not by the neglect of 
terms of high order in the eccentricity, but by the number of 
decimal places included (i.e., the truncation error is smaller 
than the round-off error). The steps in the calculation are, for 
the most part, the same as for a spherical earth; in order to 
account for the earth's ellipticity one needs, in addition, only 
to solve equation ( fTTT ) and to evaluate the series d!21 l. Since 
this approach is sufficiently convenient even for routine use, 
it is unnecessary to use an approximate method which is valid 
only for small distances. 

(The tables are shown on the following pages.) 



TABLES for computing geodesies 1. 



Arg 


log a 


-A 


log/3 


A 


log 7 A 


log a' 


-A 


log/3' 


A 


4.4 


5.314425 13 


1 


3.5124 


2000 




1 698 970 


o 


3.035 


200 


A K 
4.0 


C Q1 /I /IOC 1 O 

0.O14 4Z0 lz 


n 
U 


Q 71 O/l 
O. ( 1Z4 


oooo 
ZUUU 




T AOS 070 

i.oyo y ( u 


n 
U 


q oqp; 
o.zoO 


ooo 
zUU 


A fi 
4.0 


0.O14 4Z0 lz 


i 
1 


Q Ol O/l 

o.yiz4 


oooo 
ZUUU 




T ac\Q 070 

i.oyo y ( u 


n 
U 


q a qp; 
O.4o0 


ooo 
zUU 


A 7 
4. / 


0.O14 4Z0 11 


z 


Oil 0/1 
Z. 1 1Z4 


oooo 
ZUUU 




T ac\Q 070 

i.oyo y ( u 


n 
U 


o.DoO 


ooo 
ZUU 


A Q 
4.o 


K QI A AOK HO 
0.O14 4Z0 uy 


Q 

o 


O QI 0/1 
Z.01Z4 


oooo 
ZUUU 




T C?OQ 070 

i.oyo y ( u 


n 
U 


q qqp; 
o.ooO 


ooo 
ZUU 


4.9 


5.314 425 06 


4 


2.5124 


2000 




1.698 970 


o 


2.035 


200 


q o 
o.u 


c Qi/i /Iop: no 
0.O14 4Z0 Uz 


fi 




71 0/1 
Z. ( 1Z4 


oooo 

ZUUU 




T fiQS 070 

i .oyo y i u 


n 
U 


O OQPi 
Z.ZOO 


000 
ZUU 


Q 1 

o.l 


c qi a AO A Ofi 
0.O14 4Z4 yo 


1 o 
1U 


O Ol 0/1 

z.yiz4 


oooo 
ZUUU 




T ar\Q 070 

i.oyo y ( u 


n 
u 


o a qp; 
z.4o0 


ooo 
zUU 


Q O 

o.z 


0.O14 4Z4 00 


ID 


Til O/l 
1. 1 1Z4 


oooo 
ZUUU 




T fiOQ 070 

i.oyo y ( u 


n 
U 


o fiQp; 
z.DoO 


ooo 
ZUU 


Q Q 
O.O 


K QI /I /I O/l 7H 

O.ol4 4z4 /U 


ZO 


T QI 0/1 
1.01Z4 


oooo 
zUUU 




T fiOQ 070 

i.oyo y / u 


n 
U 


o qqp; 
z.ooO 


ooo 
zUU 


3.4 


R SI 4 424 45 


40 


1.5124 


2000 




T fiQ8 Q70 


1 


1.035 


200 


Q KC\ 
O.OU 


0.O14 4Z4 UO 





T 71 O/l 
1. / 1Z4 


ooo 
zUU 




T fiOQ OfiO 

i.oyo yoy 


n 
U 


T oqp; 
l.zoO 


OO 

zU 


Q KL1 
0.01 


c qi /i /lO/i nn 

0.O14 4Z4 UU 


fi 



T 7Q0/1 
1. ( 0Z4 


000 
ZUU 




i .oyo yoy 


n 
U 


T OPiPi 
l.ZOO 


00 
ZU 


Q PiO 
O.OZ 


C QI /I /lOQ O/l 

0.O14 4Z0 y4 





1. / 0z4 


ooo 
zUU 




T fiOQ OfiO 

i.oyo yoy 


n 
U 


T o7p; 
l.z / 


OO 

zU 


q p;q 
o.Oo 


k Qi /i /lOQ on 
O.ol4 4zo oy 





T 770/1 
1. ( ( z4 


ooo 
zUU 




T fiOQ OfiO 

i.oyo yoy 


n 
U 


T oop; 
l.zyO 


OO 

zU 


3.54 


5.314 423 83 


6 


1.7924 


200 




1 698 969 


o 


1.315 


20 


q p;p; 
o.oo 


K QI /I /lOQ 77 

O.ol4 4zo / / 


7- 
/ 


T QI O/l 
1.0lZ4 


ooo 
zUU 




T fiOQ OfiO 

i.oyo yoy 


n 
U 


T qqp; 
l.ooO 


OO 

zU 


q p;£i 
o.OD 


K QI A /IOO 70 

0.014 4Z0 / U 


7 
I 


T QQO/I 
1.00Z4 


ooo 
zUU 




T fiOQ OfiO 

i.oyo yoy 


n 
U 


l.oOO 


OO 

zU 


Q ^7 

o.o / 


c Qi A /19Q fiQ 
0.014 4Z0 DO 


i 


T QP.0/1 
1.O0Z4 


ooo 

ZUU 




i .oyo yoy 


n 
U 


T Q7Pi 
l.O/O 


00 
ZU 


q p;s 
o.Oo 


K QI /I /lOQ 

O.ol4 4zo OD 




T Q70/I 
l.O ( Z4 


ooo 
zUU 




T fiOQ OfiO 

i.oyo yoy 


n 
U 


T qop; 

i.oyo 


OO 

zU 


3.59 


5.314 423 49 


8 


1.8924 


200 




1 698 969 


o 


1.415 


20 


q fin 


K QI /I /lOQ /I 1 

0.014 4Z0 41 


Q 
O 


T oi 0/1 
i.yiZ4 


ooo 

ZUU 




i .oyo yoy 


n 
U 


T /IQPi 
1.400 


00 
ZU 


q fii 

O.Dl 


c Qi /l /lOQ QQ 
0.O14 4Z0 OO 


Q 
O 


T QQO/1 

i.yoz4 


ooo 

ZUU 




i .oyo yoy 


n 
U 


T A PiPi 
1.400 


00 
ZU 


Q fiO 
O.OZ 


c Qi A /10Q OPi 
0.014 4Z0 ZO 


Q 


T o^O/i 
i.yoz4 


ooo 
zUU 




T AOS OfiO 

i .oyo yoy 


n 
U 


T a 7p; 
1.4/0 


OO 

zU 


Q fiQ 
O.DO 


K QI /I /lOQ 1 fi 
0.O14 4Z0 ID 


1 o 
1U 


T 070/1 

i.y i Z4 


000 
ZUU 




i .oyo yoy 


n 
U 


T /IQPi 

i.4yo 


00 
ZU 


3.64 


5.314 423 06 


9 


1.9924 


200 




1 698 969 


o 


1.515 


20 


Q fip> 
O.DO 


£ QI /I /IOO 07 
0.O14 4ZZ y / 


1 1 


O Ol O/l 
U.U1Z4 


ooo 
zUU 




i .oyo yoy 


1 
1 


T PiQPi 
1.000 


00 
ZU 


O.DO 


C QI /I /lOO 

0.O14 4ZZ OD 


1 o 
1U 


O OQO/1 
U.Uoz4 


ooo 
zUU 




i .oyo yoo 


n 
U 


T PiPiPi 
1.000 


00 
ZU 


Q fi7 
O.D / 


c qi a AOO 7fi 
0.O14 4ZZ / D 


1 1 


O OP;0/1 

U.U0z4 


ooo 
zUU 




i .oyo yoo 


n 
U 


T Pi7Pi 
1.0/0 


00 
ZU 


Q fiS 
O.Oo 


C QI A AOO fix. 
0.O14 4ZZ DO 


1Z 


O 070/I 

U.U ( z4 


ooo 
zUU 




i.oyo yoo 


n 
U 


i.oyo 


OO 

zU 


3.69 


5.314 422 53 


12 


0.0924 


200 




1 698 968 


o 


1.615 


20 


Q 70 
0. / U 


£ QI /I /lOO /1 1 
0.O14 4ZZ 41 


1 Q 

lo 


Oil 0/1 
U. 1 1Z4 


ooo 
zUU 




i .oyo yoo 


n 
U 


l.DOO 


00 
ZU 


Q 71 
0. / 1 


0.O14 4zz zo 


1 A 
14 


O 1 QO/1 
U.loz4 


ooo 
zUU 




i .oyo yoo 


n 
U 


l.DOO 


00 
ZU 


Q 70 


O.ol4 4zz 14 


1 A 

14 


n i c:o/i 
U.l0z4 


onn 
zUU 




i.oyo yoo 


n 
U 


T fiTK 
l.D / 


OO 

zU 


Q 7Q 

o. (6 


O.ol4 4zz UU 


10 


n 1 TO A 

U.l f z4 


onn 
zUU 




i.oyo yoo 


n 
U 


l.Dyo 


on 
zU 


3.74 


5.314 421 85 


15 


0.1924 


200 




1 698 968 


o 


1.715 


20 


Q 7P> 
0. / 


c qi /i /ioi ^n 

0.O14 4Z1 / U 


1 fi 
ID 


O 01 O/l 
U.Z1Z4 


ooo 
zUU 




T AOS QfiQ 

i .oyo yoo 


n 
U 


T 7QPi 
1. / 00 


00 
ZU 


Q 7fi 
0. / 


C QI /I /IOI C/l 

0.O14 4Zl 04 


1 7 
1 ( 


O OQO/I 
U.zoz4 


ooo 
zUU 




i .oyo yoo 


1 
1 


T 7PiPi 
1. / 00 


00 
ZU 


Q 77 

6.(1 


K QI A /IOI Q7 

O.ol4 4zl o * 


1 Q 

lo 


n oc:0/i 
U.zoz4 


onn 
zUU 




T fiOQ Ofi7 

i.oyo yo ( 


n 
U 


T 77p; 
1. / / 


OO 

zU 


Q 7S 


c: Qi/1 /ioi in 
O.ol4 4zl ly 


1 Q 

lo 


n 070/I 
U.z r z4 


onn 
zUU 




T fiOQ Ofi7 

i.oyo yo ^ 


n 
U 


T 70p; 

i. / yo 


OO 

zU 


3.79 


5.314 421 01 


20 


0.2924 


200 




1 698 967 


o 


1.815 


20 


q so 
o.oU 


£ QI /I /IOO QI 
0.O14 4ZU ol 


OO 

zU 


O QI O/l 
U.OlZ4 


ooo 
zUU 




T aciQ o^7 

i.oyo yo ( 


n 
U 


T qqp; 
l.ooO 


OO 

zU 


Q si 
0.01 


£ QI /I /IOO fil 
0.O14 4ZU Dl 


00 

zz 


O QQO/I 
U.OOZ4 


ooo 
zUU 




T fiQS 0fi7 

i .oyo yo < 


n 
U 


1.000 


00 
ZU 


Q QO 

o.oz 


K QI A /ion QO 

O.ol4 4zU oy 


OO 

zz 


n Qco/i 
U.O0Z4 


onn 
zUU 




T ac\o ofi7 

i.oyo yo ( 


n 
U 


T Q7p; 
l.o / 


OO 

zU 


Q fiQ 
O.OO 


c qi /i /IOO 1 7 
0.O14 4ZU 1 / 


OQ 
ZO 


O Q70/I 

U.o r z4 


ooo 
zUU 




T ar\Q ofi7 

i.oyo yo ( 


n 
U 


i.oyo 


OO 

zU 


3.84 


5.314 419 94 


25 


0.3924 


200 




1 698 967 


1 


1.915 


20 


3.85 


5.314 419 69 


25 


0.4124 


200 




1.698 966 


o 


1.935 


20 


3.86 


5.314 419 44 


27 


0.4324 


200 




1.698 966 





1.955 


20 


3.87 


5.31441917 


28 


0.4524 


200 




1.698 966 





1.975 


20 


3.88 


5.314 418 89 


30 


0.4724 


200 




1.698 966 





1.995 


20 


3.89 


5.314 418 59 


31 


0.4924 


200 




1.698 966 


1 


0.015 


20 


3.90 


5.314 418 28 




0.5124 






1.698 965 




0.035 
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TABLES for computing geodesies 2. 



Arg 


logo 


-A 


log/? 


A 


log 7 


A 




logo' 


-A 


log/3' 


A 


3.90 


5.314 418 28 


32 


0.512 35 


2000 






I 


698 965 





0.035 


20 


3.91 


5.31441796 


34 


0.532 35 


2000 






1 


698 965 





0.055 


20 


3.92 


5.31441762 


35 


0.552 35 


2000 






1 


698 965 





0.075 


20 


3.93 


5.31441727 


37 


0.572 35 


2000 






1 


698 965 





0.095 


20 


3.94 


5.31441690 


39 


0.592 35 


2000 






1 


698 965 


1 


0.115 


20 


3.95 


5.31441651 


41 


0.612 35 


2000 






1 


698 964 





0.135 


20 


3.96 


5.31441610 


42 


0.632 35 


2000 






1 


698 964 





0.155 


20 


3.97 


5.314415 68 


45 


0.652 35 


2000 






1 


698 964 


1 


0.175 


20 


3.98 


5.314415 23 


47 


0.672 35 


1999 






1 


698 963 





0.195 


20 


3.99 


5.314 414 76 


48 


0.692 34 


2000 






1 


698 963 





0.215 


20 


2.00 


5.31441428 


52 


0.712 34 


2000 






1 


698 963 


1 


0.235 


20 


2.01 


5.314 413 76 


53 


0.732 34 


2000 






1 


698 962 





0.255 


20 


2.02 


5.31441323 


56 


0.752 34 


2000 






1 


698 962 





0.275 


20 


2.03 


5.314412 67 


59 


0.772 34 


2000 






1 


698 962 


1 


0.295 


20 


2.04 


5.314412 08 


61 


0.792 34 


2000 






1 


698 961 





0.315 


20 


2.05 


5.31441147 


65 


0.812 34 


2000 






I 


698 961 


1 


0.335 


20 


2.06 


5.31441082 


67 


0.832 34 


2000 






1 


698 960 





0.355 


20 


2.07 


5.31441015 


71 


0.852 34 


1999 






1 


698 960 





0.375 


20 


2.08 


5.314 409 44 


74 


0.872 33 


2000 






1 


698 960 


1 


0.395 


20 


2.09 


5.314 408 70 


77 


0.892 33 


2000 






1 


698 959 





0.415 


20 


2.10 


5.314 407 93 


81 


0.912 33 


2000 






I 


698 959 


1 


0.435 


20 


2.11 


5.314 40712 


85 


0.932 33 


2000 






I 


698 958 


1 


0.455 


20 


2.12 


5.314 406 27 


89 


0.952 33 


2000 






1 


698 957 





0.475 


20 


2.13 


5.314 405 38 


93 


0.972 33 


1999 






1 


698 957 


1 


0.495 


20 


2.14 


5.314 404 45 


98 


0.992 32 


2000 






1 


698 956 





0.515 


20 


2.15 


5.314 403 47 


102 


1.012 32 


2000 






I 


698 956 


1 


0.535 


20 


2.16 


5.314 402 45 


107 


1.032 32 


2000 






I 


698 955 


1 


0.555 


20 


2.17 


5.314 40138 


112 


1.052 32 


2000 






1 


698 954 


1 


0.575 


20 


2.18 


5.314 400 26 


117 


1.072 32 


1999 






1 


698 953 





0.595 


20 


2.19 


5.314 399 09 


123 


1.092 31 


2000 






1 


698 953 


1 


0.615 


20 


2.20 


5.314 397 86 


128 


1.11231 


2000 






I 


698 952 


1 


0.635 


20 


2.21 


5.314 396 58 


135 


1.132 31 


2000 






I 


698 951 


1 


0.655 


20 


2.22 


5.314 395 23 


141 


1.152 31 


1999 






1 


698 950 


1 


0.675 


20 


2.23 


5.314 393 82 


147 


1.172 30 


2000 






1 


698 949 


1 


0.695 


20 


2.24 


5.314 392 35 


155 


1.192 30 


2000 






1 


698 948 


1 


0.715 


20 


2.25 


5.314 390 80 


162 


1.212 30 


1999 


4.207 


40 


I 


698 947 


1 


0.735 


20 


2.26 


5.314 38918 


169 


1.232 29 


2000 


4.247 


40 


I 


698 946 


1 


0.755 


20 


2.27 


5.314 38749 


177 


1.252 29 


2000 


4.287 


40 


1 


698 945 


1 


0.775 


20 


2.28 


5.314 385 72 


186 


1.272 29 


1999 


4.327 


40 


1 


698 944 


2 


0.795 


20 


2.29 


5.314 383 86 


195 


1.292 28 


2000 


4.367 


40 


1 


698 942 


1 


0.815 


20 


2.30 


5.314 38191 


203 


1.312 28 


1999 


4.407 


40 


1 


698 941 


1 


0.835 


20 


2.31 


5.314 379 88 


213 


1.332 27 


2000 


4.447 


40 


I 


698 940 


2 


0.855 


20 


2.32 


5.314 377 75 


224 


1.352 27 


2000 


4.487 


40 


1 


698 938 


1 


0.875 


20 


2.33 


5.314 375 51 


234 


1.372 27 


1999 


4.527 


40 


1 


698 937 


2 


0.895 


20 


2.34 


5.314 37317 


244 


1.392 26 


2000 


4.567 


40 


1 


698 935 


1 


0.915 


20 


2.35 


5.314 370 73 


257 


1.412 26 


1999 


4.607 


40 


I 


698 934 


2 


0.935 


20 


2.36 


5.314 36816 


268 


1.432 25 


2000 


4.647 


40 


I 


698 932 


2 


0.955 


20 


2.37 


5.314 365 48 


281 


1.452 25 


1999 


4.687 


40 


1 


698 930 


2 


0.975 


20 


2.38 


5.314 362 67 


295 


1.472 24 


1999 


4.727 


40 


1 


698 928 


2 


0.995 


20 


2.39 


5.314 359 72 


308 


1.492 23 


2000 


4.767 


40 


1 


698 926 


2 


1.015 


20 


2.40 


5.314 356 64 




1.512 23 




4.807 




I 


698 924 




1.035 
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TABLES for computing geodesies 3. 



Arg 


logo 


-A 


log/? 


A 


log 7 


A 




logo' 


-A 


log/3' 


A 


2.40 


5.314 356 64 


323 


1.512 23 


1999 


4.807 


40 


I 


698 924 


2 


1 
1 


UoO 


20 


2.41 


5.314 353 41 


338 


1.532 22 


1999 


4.847 


40 


1 


698 922 


2 


1 
1 


Uoo 


20 


2.42 


5.314 350 03 


353 


1.552 21 


2000 


4.887 


40 


1 


698 920 


2 




075 


20 


2.43 


5.314 346 50 


371 


1.572 21 


1999 


4.927 


40 


1 


698 918 


3 


1 


095 


20 


2.44 


5.314 342 79 


388 


1.592 20 


1999 


4.967 


40 


i 


698 915 


2 


1 


115 


20 


2.45 


5.314 338 91 


406 


1.612 19 


1999 


3.007 


40 


i 


698 913 


3 


1 
1 


lor 

loo 


20 


2.46 


5.314 334 85 


425 


1.632 18 


2000 


3.047 


40 


i 


698 910 


3 


1 
1 


100 


20 


2.47 


5.314 330 60 


446 


1.652 18 


1999 


3.087 


40 


i 


698 907 


3 




175 


20 


2.48 


5.314 326 14 


466 


1.672 17 


1999 


3.127 


40 


i 


698 904 


3 


1 


195 


20 


2.49 


5.314 32148 


489 


1.692 16 


1999 


3.167 


40 


i 


698 901 


3 


1 


215 


20 


2.50 


5.314 316 59 


511 


1.712 15 


1999 


3.207 


40 


i 


698 898 


4 


1 
1 


nor 

zoo 


20 


2.51 


5.314 31148 


535 


1.732 14 


1999 


3.247 


40 


i 


698 894 


3 


1 
1 


zoo 


20 


2.52 


5.314 30613 


561 


1.752 13 


1999 


3.287 


40 


i 


698 891 


4 


1 


275 


20 


2.53 


5.314 300 52 


587 


1.772 12 


1998 


3.327 


40 


i 


698 887 


4 


1 


295 


20 


2.54 


5.314 294 65 


615 


1.792 10 


1999 


3.367 


40 


i 


698 883 


4 


1 


315 


20 


2.55 


5.314 288 50 


644 


1.812 09 


1999 


3.407 


40 


I 


698 879 


4 


1 
1 


OOO 


20 


2.56 


5.314 282 06 


674 


1.832 08 


1999 


3.447 


40 


i 


698 875 


5 


1 
1 


Q £ K 

ooo 


20 


2.57 


5.314 275 32 


705 


1.852 07 


1998 


3.487 


40 


i 


698 870 


5 


1 


375 


20 


2.58 


5.314 268 27 


739 


1.872 05 


1999 


3.527 


40 


i 


698 865 


4 


1 


395 


20 


2.59 


5.314 260 88 


774 


1.892 04 


1998 


3.567 


40 


i 


698 861 


6 


1 


415 


20 


2.60 


5.314 253 14 


810 


1.912 02 


1998 


3.607 


39 


I 


698 855 


5 


1 
1 


A Q£ 
4o0 


20 


2.61 


5.314 245 04 


848 


1.932 00 


1999 


3.646 


40 


I 


698 850 


6 


1 
1 


A KK 
400 


20 


2.62 


5.314 236 56 


889 


1.95199 


1998 


3.686 


40 


i 


698 844 


6 


1 


475 


20 


2.63 


5.314 227 67 


930 


1.971 97 


1998 


3.726 


40 


i 


698 838 


6 


1 


495 


20 


2.64 


5.314 218 37 


973 


1.99195 


1998 


3.766 


40 


i 


698 832 


6 


1 


515 


20 


2.65 


5.314 208 64 


1020 


2.01193 


1998 


3.806 


40 


I 


698 826 


7 


1 
1 


ooo 


20 


2.66 


5.314198 44 


1068 


2.03191 


1998 


3.846 


40 


I 


698 819 


7 


1 
1 


pr cr cr 

ooo 


20 


2.67 


5.314187 76 


1118 


2.05189 


1998 


3.886 


40 


i 


698 812 


8 


1 


575 


20 


2.68 


5.314176 58 


1170 


2.071 87 


1997 


3.926 


40 


i 


698 804 


7 


1 


595 


20 


2.69 


5.314164 88 


1226 


2.091 84 


1998 


3.966 


40 


i 


698 797 


9 


1 


615 


20 


2.70 


5.314152 62 


1283 


2.11182 


1997 


2.006 


40 


I 


698 788 


8 


1 
1 


ooo 


19 


2.71 


5.314139 79 


1344 


2.13179 


1998 


2.046 


40 


I 


698 780 


9 


1 
1 


Q. K A 
004 


20 


2.72 


5.314126 35 


1406 


2.15177 


1997 


2.086 


40 


i 


698 771 


9 


1 


674 


20 


2.73 


5.314112 29 


1473 


2.17174 


1997 


2.126 


40 


i 


698 762 


10 


1 


694 


20 


2.74 


5.314 097 56 


1543 


2.19171 


1997 


2.166 


40 


i 


698 752 


11 


1 


714 


20 


2.75 


5.314 082 13 


1615 


2.21168 


1997 


2.206 


40 


I 


698 741 


10 


1 
1 


1 o4 


20 


2.76 


5.314 065 98 


1690 


2.23165 


1996 


2.246 


40 


I 


698 731 


12 


1 
1 


(04 


20 


2.77 


5.314 049 08 


1771 


2.25161 


1997 


2.286 


40 


i 


698 719 


11 


1 


774 


20 


2.78 


5.314 03137 


1853 


2.271 58 


1996 


2.326 


40 


i 


698 708 


13 


1 


794 


20 


2.79 


5.314 012 84 


1941 


2.29154 


1996 


2.366 


39 


i 


698 695 


13 


1 


814 


20 


2.800 


5.313 993 43 


1004 


2.31150 


998 


2.405 


20 


i 


698 682 


6 


-i 
1 


QQ A 

oo4 


10 


2.805 


5.313 983 39 


1028 


2.32148 


998 


2.425 


20 


I 


698 676 


7 


-i 
1 


Q A A 
544 


10 


2.810 


5.313 97311 


1051 


2.33146 


998 


2.445 


20 


i 


698 669 


7 


1 


854 


10 


2.815 


5.313 962 60 


1076 


2.34144 


998 


2.465 


20 


i 


698 662 


7 


1 


864 


10 


2.820 


5.313 95184 


1101 


2.35142 


998 


2.485 


20 


i 


698 655 


8 


1 


874 


10 


2.825 


5.313 940 83 


1127 


2.36140 


997 


2.505 


20 


I 


698 647 


7 


1 


884 


10 


2.830 


5.313 929 56 


1152 


2.371 37 


998 


2.525 


20 


I 


698 640 


8 


1 


894 


10 


2.835 


5.313 918 04 


1180 


2.38135 


998 


2.545 


20 


i 


698 632 


8 


1 


904 


10 


2.840 


5.313 906 24 


1207 


2.39133 


997 


2.565 


20 


i 


698 624 


8 


1 


914 


10 


2.845 


5.313 89417 


1234 


2.401 30 


998 


2.585 


20 


i 


698 616 


8 


1 


924 


10 


2.850 


5.313 88183 




2.41128 




2.605 




I 


698 608 




1 


934 
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TABLES for computing geodesies 4. 



Arg 


log a 


—A 


log/? 


A 


log 7 


A 


log a' 


-A 


log/?' 


A 


2.850 


5.313 88183 


1264 


2.411279 


9974 


2.605 


20 


1.698 608 


8 


1.934 


10 


2.855 


5.313 869 19 


1293 


2.421253 


9974 


2.625 


20 


1.698 600 


9 


1.944 


10 


2.860 


5.313 856 26 


1323 


2.431227 


9974 


2.645 


20 


1.698 591 


9 


1.954 


10 


2.865 


5.313 843 03 


1353 


2.441201 


9973 


2.665 


20 


1.698 582 


9 


1.964 


10 


2.870 


5.313 829 50 


1385 


2.451 174 


9972 


2.685 


20 


1.698 573 


9 


1.974 


10 


2.875 


5.313 815 65 


1417 


2.461 146 


9972 


2.705 


20 


1.698 564 


10 


1.984 


10 


2.880 


5.313 80148 


1450 


2.471 118 


9971 


2.725 


20 


1.698 554 


9 


1.994 


10 


2.885 


5.313 786 98 


1484 


2.481089 


9970 


2.745 


20 


1.698 545 


10 


2.004 


10 


2.890 


5.313 772 14 


1518 


2.491059 


9970 


2.765 


20 


1.698 535 


10 


2.014 


9 


2.895 


5.313 756 96 


1553 


2.501029 


9969 


2.785 


19 


1.698 525 


11 


2.023 


10 


2.900 


5.313 74143 


1590 


2.510 998 


9968 


2.804 


20 


1.698 514 


10 


2.033 


10 


2.905 


5.313 725 53 


1626 


2.520 966 


9968 


2.824 


20 


1.698 504 


11 


2.043 


10 


2.910 


5.313 709 27 


1664 


2.530 934 


9966 


2.844 


20 


1.698 493 


11 


2.053 


10 


2.915 


5.313 692 63 


1702 


2.540 900 


9966 


2.864 


20 


1.698 482 


11 


2.063 


10 


2.920 


5.313 675 61 


1742 


2.550 866 


9965 


2.884 


20 


1.698 471 


12 


2.073 


10 


2.925 


5.313 658 19 


1783 


2.560 831 


9965 


2.904 


20 


1.698 459 


12 


2.083 


10 


2.930 


5.313 640 36 


1824 


2.570 796 


9963 


2.924 


20 


1.698 447 


12 


2.093 


10 


2.935 


5.313 622 12 


1866 


2.580 759 


9963 


2.944 


20 


1.698 435 


12 


2.103 


10 


2.940 


5.313 603 46 


1909 


2.590 722 


9962 


2.964 


20 


1.698 423 


13 


2.113 


10 


2.945 


5.313 584 37 


1953 


2.600 684 


9961 


2.984 


20 


1.698 410 


13 


2.123 


10 


2.950 


5.313 564 84 


1999 


2.610 645 


9960 


1.004 


20 


1.698 397 


13 


2.133 


10 


2.955 


5.313 544 85 


2045 


2.620 605 


9959 


1.024 


20 


1.698 384 


14 


2.143 


10 


2.960 


5.313 524 40 


2093 


2.630 564 


9958 


1.044 


20 


1.698 370 


14 


2.153 


10 


2.965 


5.313 503 47 


2141 


2.640 522 


9957 


1.064 


19 


1.698 356 


14 


2.163 


10 


2.970 


5.313 482 06 


2191 


2.650 479 


9956 


1.083 


20 


1.698 342 


15 


2.173 


10 


2.975 


5.313 460 15 


2241 


2.660 435 


9956 


1.103 


20 


1.698 327 


15 


2.183 


10 


2.980 


5.313 437 74 


2293 


2.670 391 


9954 


1.123 


20 


1.698 312 


15 


2.193 


10 


2.985 


5.31341481 


2347 


2.680 345 


9953 


1.143 


20 


1.698 297 


16 


2.203 


9 


2.990 


5.313 39134 


2400 


2.690 298 


9952 


1.163 


20 


1.698 281 


15 


2.212 


10 


2.995 


5.313 367 34 


2457 


2.700 250 


9951 


1.183 


20 


1.698 266 


17 


2.222 


10 


1.000 


5.313 342 77 


2513 


2.710 201 


9950 


1.203 


20 


1.698 249 


17 


2.232 


10 


1.005 


5.31331764 


2571 


2.720 151 


9948 


1.223 


20 


1.698 232 


17 


2.242 


10 


1.010 


5.313 29193 


2631 


2.730 099 


9948 


1.243 


20 


1.698 215 


17 


2.252 


10 


1.015 


5.313 265 62 


2691 


2.740 047 


9946 


1.263 


19 


1.698 198 


18 


2.262 


10 


1.020 


5.313 238 71 


2754 


2.749 993 


9945 


1.282 


20 


1.698 180 


18 


2.272 


10 


1.025 


5.313 211 17 


2818 


2.759 938 


9943 


1.302 


20 


1.698 162 


19 


2.282 


10 


1.030 


5.313 182 99 


2883 


2.769 881 


9943 


1.322 


20 


1.698 143 


19 


2.292 


10 


1.035 


5.313 15416 


2949 


2.779 824 


9941 


1.342 


20 


1.698 124 


20 


2.302 


10 


1.040 


5.313 124 67 


3018 


2.789 765 


9939 


1.362 


20 


1.698 104 


20 


2.312 


10 


1.045 


5.313 094 49 


3087 


2.799 704 


9939 


1.382 


20 


1.698 084 


20 


2.322 


10 


1.050 


5.313 063 62 


3159 


2.809 643 


9936 


1.402 


20 


1.698 064 


21 


2.332 


10 


1.055 


5.313 032 03 


3232 


2.819 579 


9936 


1.422 


20 


1.698 043 


22 


2.342 


9 


1.060 


5.312 999 71 


3306 


2.829 515 


9934 


1.442 


19 


1.698 021 


22 


2.351 


10 


1.065 


5.312 966 65 


3383 


2.839 449 


9932 


1.461 


20 


1.697999 


22 


2.361 


10 


1.070 


5.312 932 82 


3460 


2.849 381 


9931 


1.481 


20 


1.697977 


23 


2.371 


10 


1.075 


5.312 898 22 


3541 


2.859 312 


9929 


1.501 


20 


1.697 954 


24 


2.381 


10 


1.080 


5.312 862 81 


3623 


2.869 241 


9928 


1.521 


20 


1.697 930 


24 


2.391 


10 


1.085 


5.312 826 58 


3706 


2.879 169 


9926 


1.541 


20 


1.697906 


25 


2.401 


10 


1.090 


5.312 789 52 


3791 


2.889 095 


9924 


1.561 


20 


1.697881 


25 


2.411 


10 


1.095 


5.312 75161 


3879 


2.899 019 


9922 


1.581 


19 


1.697856 


26 


2.421 


10 


1.100 


5.312 712 82 




2.908 941 




1.600 




1.697 830 




2.431 





